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Abstract 

The dynamics of the highly nonlinear fifth order KdV-type equation is dis- 
cussed in the framework of the Lagrangian and Hamiltonian formalisms. The 
symmetries of the Lagrangian produce three commuting conserved quantities 
that are found to be recursively related to one-another for a certain specific 
value of the power of nonlinearity. The above cited recursion relations are 
obeyed with a second Poisson bracket which sheds light on the integrability 
properties of the above nonlinear equation. It is shown that a Miura-type 
transformation can be made to obtain the fifth order mKdV-type equation 
from the fifth order KdV-type equation. The spatial dependence of the fields 
involved is, however, not physically interesting from the point of view of the 
solitonic solutions. As a consequence, it seems that the fifth order KdV- and 
mKdV-type equations are completely independent nonlinear evolution equa- 
tions in their own right. 
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1 Introduction 



In the classical soliton theory, the ideas of integrability, elastic scattering of soli- 
tons, stability and preservation of their shape, etc., are fairly well established to be 
very intimately connected. In recent years, the stability of the soliton solutions to 
the highly nonlinear fifth order KdV-type equations has been studied extensively 
[1-5]. In particular, the effect of the dispersion due to higher order derivative terms 
on solitons has been studied for the fifth order KdV- and mKdV-type equations 
which are useful in the context of fluid mechanics and plasma physics [2-4]. It has 
been shown that the nonzero dispersion effects lead to the narrowness of solitonic 
solutions [5]. Furthermore, it has been established that the solitons are stable for 
p < 8 where p corresponds to the power of nonlinearity [6] (see, e.g., eqn. (1) below). 
The analytic expression for the soliton solution to the fifth order KdV(FOKdV)- 
type equation and two analytic solutions to the fifth order mKdV(FOmKdV)-type 
equations have also been found [7]. 

It is obvious that the above equations, their stability and solutions are of consid- 
erable theoretical interests for the understanding of some physical processes in the 
realm of fluid mechanics and plasma physics. 

It is a well known fact that 1+1 dimensional (super) integrable equations (e.g., 
(super)KdV and (super) Boussinesq, etc.) play a prominent role in the understand- 
ing of (super)string theory, W-type (super) algebras, statistical mechanics at critical 
temperature, 2D gravity theories, fluid mechanics, etc. Their integrability proper- 
ties stem from the presence of an infinite number of commuting conserved quantities. 
The existence of second Hamiltonian structure, zero- curvature representation, Lax- 
pair formulation, multi-soliton solutions, Miura-type transformations, etc., are also 
some of the key features of integrability (see, e.g., [8-10]). The purpose of the 
present work is to study the highly nonlinear FOKdV and FOmKdV equations 
in the light of the above integrability properties as it is always worthwhile to look 
for a nonlinear evolution equation that is integrable or partially integrable [11]. To 
this end in mind, we first develop a Lagrangian formulation for the above dynam- 
ical equations following the procedure adopted by Cooper etal.[12] in the context 
of "compactons" and derive conserved quantities from the symmetries of the La- 
grangian due to Noether theorem. We discuss a consistent Hamiltonian formulation 
for the above systems and concentrate on the Poisson bracket structures to study 
the commuting properties of the conserved quantities. It turns out that for a specific 
value of the power of nonlinearity, the conserved quantities are recursively related 
to one-another. However, the second Poisson bracket structure is found to be not 
a consistent bracket for an integrable system [13]. The FOmKdV- type equation is 



2 



derived from the FOKdV-type equation by a Miura-type transformation. It turns 
out that the spatial dependence of the fields involved is such that the meaning of soli- 
tonic solutions to these equations is lost. In contrast to the usual KdV and mKdV 
equations which are related by Miura maps, this result implies that FOKdV and 
FOmKdV equations are independent in their own right. A consistent Lagrangian 
and Hamiltonian formulations for the FOmKdV equation is also developed. 

The material of our work is organized as follows. In Sec. 2, we derive the FOKdV 
equation from the scalar field Lagrangian and discuss its symmetry properties. Sec. 
3 is devoted to the Hamiltonian structures and discussion of the integrability proper- 
ties. We try to understand how far and how much these properties can be discussed 
vis-a-vis the usual KdV and mKdV equations. In Sec. 4, we derive FOmKdV 
equation from the starting FOKdV equation by a Miura-type transformation and 
demonstrate that such type of relation entails the spatial dependence of the fields 
involved to be some fractional power over spatial variable(x). As a consequence, it 
is difficult to provide some physical interpretation to this result in the language of 
solitonic solution. Finally, in Sec. 5, we make some concluding remarks and provide 
an outlook of the whole discussion. 



2 Lagrangian Formulation 

We derive here the FOKdV-type equation 

ut + 1 u 5x + (3 u 3x + a u v u x = 0, (1) 

from the first-order Lagrangian density C that is expressed in terms of the derivatives 
(0t = dt<p, 4> x = d x <fi) on a scalar field <p. The continuous symmetries of the action 
(S — J dt J dx C) lead to the derivation of three conserved quantities due to Noether 
theorem. We begin with the following first-order Lagrangian 

l ± ± , a (<p x y +2 (5 l2 , 7 
(p+l)(p + 2) 

where p > and c-number parameters a, (3, 7 of (1) and (2) obey certain restrictions 
for the physically meaningful solutions to the equation (1) in the context of fluid 
mechanics and plasma physics (see, e.g., [4]). The Euler-Lagrange equation 

d 5L d 5L 

with 

SL 



l w = / ix c . / a* ^ <j> t + (p a + ( ?f v + * 2) - 1 A + 1 A] • (2) 



TT - Vx 

n - y , 



8<f> x 2 (p + 1) 
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leads to the following dynamical equation of motion 



>:,:t 



+ a (<p x ) p 4> 2x + (3 <p4x + 7 <Pgx = 0. (5) 



The identification 4> x = u(x,t) leads to the derivation of (1). 

We would like to dwell a bit more on the symmetries of the above Lagrangian as 
they provide a physical insight into the system described by it. Since only derivatives 
are present in the above Lagrangian, it is straightforward to check that the following 
transformation with a global infinitesimal parameter e 

M = 2e, (6) 

is a symmetry transformation because 5\C = 0. Furthermore, the space and time 
translations x ^ x — £, t — > i — e lead to the following transformations on the scalar 
field (j) ( since <p(x', f) = <p(x, t) ) 

S 2 (j) = e (j) x , 

5 3 <f) = e<f) t . (7) 

These are also the symmetry transformations of the action S as the above Lagrangian 
undergoes the following change 

d 

5 2 C = eC x = —(eC), 

hC = eC t = ^{eC), (8) 
where e is a global infinitesimal parameter (as was the case in (6)). 

The continuous symmetry transformations (6) and (7) imply presence of con- 
served quantities (I's) due to well-known Noether theorem. These are juxtaposed, 
with the identification (fi x = u, as follows 

Ii — I dxcj) x — > / dxu, 
/3 = J dx k^-(v+l)(v + 2)-2^^ 



J dx(^ 



(p+l)(p + 2) 

)• (9) 



(5 2 a u p+2 7 2 

2 Ux ~ (p+l)(p + 2) " 2 U2x 



Note that in the derivation of ^3, we have used 

'5 3 (j) dC 



1, = /**, = C), (10) 
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as the Lagrangian transforms to the time derivative of itself under global time trans- 
lation (7). It can be easily seen that the time translation generator "energy" (1%) 
is nothing but the Hamiltonian function (H(p)) corresponding to the first-order 
Lagrangian function (2). It becomes transparent in the following Legendre transfor- 
mation 

Hip) = [ dx[U^ 4> t \ - L{p) = h. (11) 



The conservation laws for the quantities lip, corresponding to the "area" under 
u and "momentum", can be understood directly from the form of the dynamical 
equation. It is obvious from (1) that I\ is a conserved quantity due to the fact that 
r.h.s. of this equation is a total space derivative. For the conservation law of the 
space translation generator (J 2 ), the equation of motion (1), implies that 



uu t = —z- 

ox 



d r ct u p+2 1 1 i 

+ P (uu 2x - -u 2 x ) + 7 (uu 4x - u x u 3x + -u\ x ) . (12) 



(p + 2) ' " ^ 2' 



In general, for an equation of the type (1) with any order of odd derivatives on u 
and a nonlinear term of the type u p u x , it can be seen that u Ut is always a total 
space derivative. It is because of the fact that for n — 0, 1, 2, 3 

d n ~ 1 1 

U(2n+l)xU=—[ J2 ( _1 ) r U (2n-r)x U rx + (- 1)™- (w n:E ) 2 ] (13) 
OX r=0,n>r 1 

where u 0x — u, r — 0, 1, 2, 3, (n — 1). This demonstrates that the "momen- 
tum" I 2 is always a conserved quantity for the kind of equation we are discussing. 

From the above Hamiltonian function, one can derive the dynamical equation (5) 
(and (1) with proper identification) by exploiting the following canonical Poisson- 
bracket structure for the scalar field 

{<Pix,t),<P y iy,t)}% =6(x-y), (14) 

in the Hamilton equation of motion 

<k = {0, H{p)}% ee - 7 5x - p <p 3x - a (15) 

It is evident that the total space derivative of the above equation is nothing but 
the dynamical equation (5) which, in turn, takes the form of equation (1) with the 
identification (fi x = u. Similar operation on the Poisson bracket (14) leads to the 
Hamiltonian structure on u (the first PB) as 

{«(*), «(!/)}« =^S(x-y). (16) 

Thus, we establish the fact that the Lagrangian and Hamiltonian formulations for 
the scalar field lead to the derivation of the Poisson-bracket structure for the u fields 
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as has been taken by Karpman in Ref. [4] for the Hamiltonian function J 3 . It is 

interesting to point out that all the conserved quantities of equation (9) commute 
with one-another under the choice of the Poisson bracket (16). 

3 Second-Hamiltonian Structure and Recursion Re- 
lations 



Here we demonstrate that for p = 1, these conserved quantities are recursively 
related to one-another for any arbitrary values of a, (3, 7 as there exists a second 
Poisson bracket for the equation of motion (1) with the Hamiltonian as the "mo- 
mentum" expression I 2 . Thus, to some extent, the FOKdV-type equation does 
mimic the usual KdV equation. 



We know that the KdV equation is an integrable equation because it supports an 
infinite number of conserved quantities that are recursively related to one-another 
due to the existence of a "consistent" second Hamiltonian structure [13]. The first- 
and the second brackets satisfy the famous Jacobi identities. Here, we try to see how 
much and how far, the FOKdV equation can mimic the anologue of the recursion 
relations of KdV equation. To this end, it can be seen that the following Hamilton 
equations 

u t = {u,h}V={u,hY 2 \ (17) 
with the first bracket (cf., equation (16)) 

{u(x,t),u(y,t)}^ = -D^5(x-y), (18) 

and the second bracket 

{u(x,t),u(y,t)}W = D^S(x-y), (19) 



where, 



dy 

D (2) = [ 7 d 5 + /3<9 3 + Cu p 3 + {(p+l) C-ayu^Uy], (20) 

lead to the derivation of equation (1) for any arbitrary value of the constant C. 
For the FOKdV-type equation, the analogue of KdV type recursion relations, with 
Poisson brackets (18) and (19), are 

DW(^)=i?G0(%i), (21) 

where k > 2 and opertators D^ 1 ' 2 ) can be computed from the brackets (18) and (19). 
These operators for our discussion are 

D( 2 ) = -[ 7 d 5 + /3<9 3 + Cw p <9- (C-a)u p ~ l du], 

D« = d. (22) 
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For k = 2, the above recursion relation leads to the following expression for I 2 

However, the comparison with the starting momentum expression ( i.e. I 2 = 
I / dx u 2 of equation (9)), puts certain restrictions on p and C as given below: 

p=l and C = l + a. (24) 

As a result, the operator ZK 1 ) remains intact but the operator £)( 2 ) modifies as 
follows: 

D \2) = + p& + (\ + a )ud-du ]. (25) 

It will be noticed that the form of the equation (1) is changed because there is a 
restriction on p ( i.e. p — 1 ) but c-number parameters a, (3 and 7 remain intact as 
there are no restriction on them due to the recursion relations. 

We go a step further to compute I 3 from the above recursion relation with the 
modified version of second Poisson bracket (19). The corresponding recursion rela- 
tion with the above DP) (cf. (25)) 

leads to 

/3 = /dx(^-|« 3 -|«L). (27) 

The comaprison with the "energy" expression I 3 (for p = 1 in equation (9)), puts no 
constraint on a, (3, 7 and all the three conserved quantities are recursively related. 
It can be checked that, for p = 1, all of them commute with one-another under the 
Poisson bracket structure (19) as well. 

The next order recursion relation should yield the next conserved quantity I4 if 
the second Hamiltonian structure is a "consistent" bracket of an integrable equation. 
In fact, the following expression comes out from the next recursion relation for any 
arbitrary value of a, (5 and 7: 

SI a <y 2 

-T~ = l( u lx + uu 4x - 2u x u 3x ) + -u 3 + (3(u u 2x - u 2 x ) + P 2 u^ x + — u 3 

OU 6 

1 7 

+ l 2 u 8x + 2Pl u 6x + ap(-u 2 x + 2uu 2x ) + a-f(-u 2 2x + 2uu 4x + 3u x u 3x )(28) 

To obtain the expression for I4, one should see that the r.h.s. is a total variation of 
some quantity. It turns out that for the choice a — 3, the r.h.s. reduces to a simpler 
form as given below 

r 23 1 
SI 4 = [7^(—ul x + uu 4x + u x u 3x ) + 7P(uu 2x + -u 2 x ) 

7 , 

+ P 2 u ix + -f 2 u 8x + 2« 3 + 2 /^7 uq x \ Su. (29) 
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Except the first term, all the other terms are total variation on a certain local 
quantity modulo some total space derivatives. This can be seen as follows: 

11 1 

(-2 U l + u u 2x)Su = «2x) u 8x 6u = -S(ul x ) 

1 17 7 

u 6x 5u = --6{u 2 3x ), u 4x Su = -5(ul x ), -u 3 5u = - 5(u 4 ). (30) 

It is interesting to note that the first term can not be cast as a total variation of a 
definite local quantity. The argument runs as follows. There are only four candidates 
with three u fields and four derivatives on them 

u 2 u 4x ,ul x u,u u x u 3x , v 2 u 2x , (31) 

whose variation and/or the variation of their linear combination might produce the 
combination of the first term. It turns out that the variation on the last term is zero 
modulo some total space derivative term. All the other three candidates produce 
the same result due to the variation on them. For instance, for one of them the 
variation is 

(2u u Ax + 4u x u 3x + 3ul x )5u = - 5(u 2 u 4x ). (32) 

Thus, the FOKdV equation is not an integrable equation because the second Hamil- 
tonian structure is not consistent [13]. It can be also seen that the second- Poisson 
bracket does not satisfy the Jacobi identities (see, e.g., for more detail, Ref. [12] for 
a similar kind of situation). 

4 Towards a Miura-type transformation 



In this Sec, we derive the FOmKdV-type equation (i.e., p = 2) from the starting 
FOKdV equation by Miura-type transformation (see, e.g., [13]). This equation is 
also of theoretical interests in the context of fluid mechanics and plasma physics. 
We develop here the Lagrangian formulation for the FOmKdV and discuss in a 
nut-shell the Poisson bracket structures for this system on similar lines as that of 
the FOKdV equation ( cf. Sec. 2). The solitonic solution to this equation has been 
found in Ref. [7] where two analytic solutions have been obtained. 

To start with, we know that for the usual KdV equation, there exists a trans- 
formation that relates the original KdV equation to the mKdV. The question 
we address is : Can one derive the fifth-order mKdV equation from the starting 
equation (1)? To answer this question we make the following transformation 

u(x,t) = v x (x,t) - — v 2 (x,t), (33) 

in equation (1). It will be noticed that the inclusion of higher order derivative 
terms on the r.h.s. of equation (33) leads to more complicated expressions which 



8 



are difficult to handle. The ensuing equation is realized on v(x, t) as 

a 2 

v t + 7 v 5x + (3 v 3x - — v 2 v x = 0, (34) 

which can be recognized as the FOmKdV and can be re-expressed in the form of a 
conservation law as given below: 



d 
ox 



a 



iv 4x + fi v 2x - — v d \ = 0. (35) 



It should be added, however, that for the validity of (35), the field v(x,t) has to 
satisfy an extra condition given by 

v x v ix + 2 v 2x v 3x = 0. (36) 

We have to find a solution for the above restriction if Miura transformation (33) has 
to make some sense. The equation (36) can be re-written as a total space derivative 
of a certain quantity, given by: 

-^[v x v 3x + ^(v 2x ) 2 ] = 0. (37) 

The quantity in the square bracket has to be a constant. Choosing the constant to 
be zero, we obtain the x dependence of v as 

v(x) = A x 5/3 , (38) 

for an arbitrary constant A. It will be noticed that equation (38) is a nontrivial 
solution for the Miura transformation and there is no constraint on the t-dependence. 
It is also evident that for such a choice (v(x) ~ x 5 / 3 ), the spatial dependence of u 
field will also change accordingly. Thus, there is a problem about solitonic solutions 
for such kind of restriction. It is well known that for the usual KdV equation, no 
such restrictions as (36-38) arise. For instance, if 7 = and p — 1 in the starting 
equation (1), there will be no restrictions like (36-38) and for any arbitrary a, (3, we 
have 

u t + (3 u 3x + a uu x = 0, 
u(x,t) = v x (x,t) - — v 2 , 

v t + /3v 3x -^ v 2 v x = 0, (39) 
as KdV equation, Miura transformation and the mKdV equation. 

We conclude from equation (38) that FOKdV and FOmKdV equations are in- 
dependent in their own right. There is no Miura type of mapping which can relate 
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these equations preserving the nature of stationary solitonic solutions that have been 
obtained in Ref. [7]. 



The FOmKdV equation of motion (34) can be derived from the following first- 
order Lagrangian 



dx 



\<t>* h - ^(^) 4 - § (M 2 + \ (M 2 )- (40) 



The Euler-Lagrange equation, emerging from the invariance of the action (SS = 0), 
leads to the following dynamical equation: 



" 2 , [2 



(f>xt ~ -^{4>x) <p2x + P <p4x + 7 06x = 0. (41) 
O jJ 

The identification <p x = v(x, t) leads to the rederivation of (34). The anologue of the 
symmetry transformations (6) and (7) for the above Lagrangian lead to the following 
conserved currents due to the Noether theorem 



Ji = J dx v, J2 = - J dx v 2 , 

j 3 - iH^+fr-H)- (42) 



The equation of motion (34) can be derived from the Hamilton equation as well. The 
corresponding Hamiltonian functions are none other than J 3 and J 2 as illustrated 
below: 

vt = {v,Js} w ={v,J 2 y 2) , (43) 
where the first bracket and the second brackets are 

{v(x,t),v(y,t)}V = -0- 5 (x-y), 

2 

{v(x,t),v(v,t)} {2) = hd 5 + pd 3 + Cv 2 d+(3C + -^)v y v}S(x-y), (44) 
for arbitrary constant C and d — J^. 

5 Discussions 



The integrability properties play a pivotal role in furnishing the wealth of infor- 
mations on the rich mathematical structure encoded in a given nonlinear evolution 
equation. This, in turn, leads to gain some deep physical insight into the phenom- 
ena that are described by these equations. In the light of this argument, we have 
studied the FOKdV and FOmKdV equations. We have developed a Lagrangian 
formulation for the FOKdV equation and obtained the conserved quantities from 
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the symmetry properties of this Lagrangian. We have shown that these conserved 
quantities are recursively related to one-another for p = 1 due to the existence of 
bi-Hamiltonian structures. The recursion relations do not lead to the determination 
of more than three conserved quantities. Thus, as it appears, the FOKdV equation 
is not a completely integrable system. 

The FOKdV equation mimics some of the key features of the partially integrable 
system. For instance, the soliton solutions are stable and preserve their shape for 
considerably large values of the power of nonlinearity. There exists a Miura type 
transformation that relates FOKdV and FOmKdV equations for the specific space- 
dependence of the fields on which these equations are realized. However, from the 
point of view of the solitonic solutions, this dependence is not an interesting feature. 
This leads to the conclusion that FOKdV and FOmKdV equations are independent 
in their own right and there exists no Miura-type mapping which can relate these 
equations without spoiling the solitonic solutions to these equations. Furthermore, 
for jo = 1, the conserved quantities are recursively related to one-another (due to 
the existence of two Poisson bracket structures). One can also see that the FOKdV 
and FOmKdV equations can be expressed as an Abelian-type zero- curvature rep- 
resentation (F tx = 0) 

F tx = d t A x - d x A t = [- + A t , — + A x \ = 0, (45) 

where gauge connection A x and A t can be read off from the FOKdV and FOmKdV 
equations respectively, as listed below 



A x = u(x,t) and v(x,t), 

u p+1 a 2 

A t = -(7M 4 x + Pu 2x + -) and — - v 3 - f3v 2x - -yv 4x . (46) 

p + 1 lop 

It will be an interesting venture to study the integrability properties of these 
equations where we can use the Painleve test [11] and consider the multi-soliton 
solution by exploiting the Hirota's bilinear form and its generalization [14]. It is 
worthwhile to mention that in Ref. [14], it has been claimed that the existence of 
three solitonic solutions imply integrability. These are the issues for future investi- 
gations and further works have to be done in this direction. 
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